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Abstract. We review our recent study of the pion mass dependence of the p and a resonances, 
generated from one-loop SI/ (2) Chiral Perturbation Theory (ChPT) with the Inverse Amplitude 
Method (I AM). In order to properly account for the Adler zero region, we also review the recently 
obtained modified version of the IAM; which is based on analyticity, elastic unitarity and ChPT at 
low energies, thus yielding the correct pion mass dependence of the resonance pole positions up to 
next-to-leading order in ChPT. As main results we find that the p %% coupling constant is almost 
m K independent and that M p shows a smooth m K dependence while that of the a shows a strong 
non-analyticity. These findings are important for studies of the meson spectrum on the lattice. 
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INTRODUCTION 

Light hadron spectroscopy lies beyond the realm of perturbative QCD. Although lat- 
tice QCD provides, in principle, a rigorous way to extract non-perturbative quanti- 
ties, present calculations use relatively high quark masses, at least for studies of scalar 
mesons, i. e., OL0]- Thus, appropriate extrapolation formulas are called for. Chiral Per- 
turbation Theory (ChPT) [3] provides such extrapolations, since it is built as an expan- 
sion in momenta and masses, generically 0(p/4nf n ) 2 , of a Lagrangian involving the 
Goldstone Bosons of the QCD chiral symmetry breaking (pions), compatible with all 
QCD symmetries. ChPT is renormalized order by order by absorbing loop divergences 
in the parameters of higher order counterterms (low energy constants - LECs), which 
are the coefficients of the energy and mass expansion, so they have no quark mass de- 
pendence. Their values depend on the QCD dynamics, and have to be determined from 
experiment. In SU(2)-ChPT %% scattering only four LECs appear, denoted 

The ChPT expansion provides a systematic and model independent description of 
how the observables depend on QCD parameters, like the quark masses, and this can be 
implemented systematically up to the desired order in the ChPT expansion. 

We review here our recent derivation of a modified version of the IAM ; based on 
dispersion theory, unitarity and ChPT to next-to-leading order (NLO), which we use to 
predict the quark mass dependence of the a and p mesons 

We focus only on the o and the p, so it is enough to work with the lightest quarks u, d 
in the isospin limit with a mass m = (m u + m^)/2. Since m n is given by m\ ~ m + ... 
[|3|], studying the m dependence is equivalent to study the m % dependence. 



UNITARIZATION AND DISPERSION THEORY 



The o and p resonances appear as poles in the second Riemann sheet of the (/, J) = (0, 0) 
and (1,1) partial waves of isospin / and angular momentum J respectively. For these 
partial waves, elastic unitarity implies, for physical values of s: 

Imt(s) = o(s)\t(s)\ 2 ^Im (t(s)- 1 ) =-o(s), with o(s) = 2p/^fs, (1) 

where s is the Mandelstamm variable and p is the center of mass momentum. Conse- 
quently, the imaginary part of the inverse amplitude is known exactly. ChPT amplitudes, 
being an expansion t = t 2 + t 4 H with = 0(p k ), satisfy Eq. (OQ) just perturbatively: 

lmt 4 (s) = o(s)\t 2 (s)\ 2 1 => lmt 4 (s)/t 2 (s) 2 = o(s) 1 (2) 

and cannot generate poles. Therefore the resonance region lies beyond the reach of 
standard ChPT. However, it can be reached by combining ChPT with dispersion theory 
either for the amplitude [6] or the inverse amplitude through the IAM 0,[a@]. 

The elastic IAM [4] uses the ChPT series and elastic unitarity to evaluate a dispersion 
relation for the inverse amplitude. The analytic structure of l/t consists on a right cut 
from threshold to °°, a left cut from — °° to s = 0, and possible poles coming from zeros 
of t. For scalar waves, t vanishes at the so called Adler zero, sa, that lies on the real 
axis below threshold, thus within the ChPT region of applicability. Its position can be 

obtained from the ChPT series, i.e., sa = s 2 + s 4 H , where t 2 vanishes at s 2 , h + £4 at 

s 2 + s 4 , etc. 

We write then a once subtracted dispersion relation for the inverse amplitude, where 
we have chosen the subtraction point to be the Adler zero: 

t(s) % Jrc (s' -sa){s'-s) 

where "LC" stands for a similar integral over the left cut and "PC" stands for the 
contribution of the pole at the Adler zero. Since t 2 is real on the real axis and £4 has 
the same analytic structure as £,we can similarly write 

i 4« to, t WMtf (tt/ti )+PC (t4/2) (4) 

t 2 {s) 2 % Jrc (s'-s 2 )(s'-s) 

where we have now subtracted at s 2 , which is the LO approximation to the Adler 
zero, and PC stands for the contribution of the pole at s 2 . We can now use unitarity, 
Eqs. (OQ) and ©, to find that the imaginary parts on the right cut of both dispersion 
relations are exactly opposite to each other. Since the LC integral is weighted at low 
energies, we can use ChPT to approximate LC( l/t) ~ — LC{t 4 /t 2 )- The pole contribution 
PC (l/t) can also be evaluated with ChPT since it involves derivatives of t evaluated at 
the Adler zero, where ChPT is perfectly justified. Finally, we approximate with ChPT 
(s — sa)(s' — sa) — (s - s 2 )/(s' - s 2 ). Altogether, we find a modified IAM (mlAM) 
formula: 

mlAM = l l A mIAM = , v {s 2 -S A ){s-S 2 )[t' 2 (s 2 )-t' A {s 2 )] 

t 2 -t 4 +A mIAM1 41 21 s-s A 
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FIGURE 1. Left: Movement of the a (dashed lines) and p (dotted lines) poles for increasing m K 
(direction indicated by the arrows) on the second sheet. The filled (open) boxes denote the pole positions 
for the a (p) at pion masses m K = 1,2, and 3 x m^ hys , respectively. For m K — 3m£ hys three poles 
accumulate in the plot very near the threshold. Note that all poles are always far from the Adler 
zero (circles). Right: The grey band shows the m K dependence of p pole mass from the IAM versus 
recent lattice results from [1]. The dashed line is the IAM result for N c = 10. 



The standard IAM is recovered for A = 0, which holds exactly for all partial waves 
except the scalar ones. In the original IAM derivation [0, Hi A mIAM was neglected, since 
it formally yields a NNLO contribution and is numerically very small, except near the 
Adler zero, where it diverges. However, if A mlAM is neglected, the IAM Adler zero 
occurs at S2, correct only to LO, it is a double zero instead of a simple one, and a spurious 
pole of the amplitude appears close to the Adler zero. All of these caveats are removed 
with the mlAM, Eq. ©. The differences in the physical and resonance region between 
the IAM and the mlAM are less than 1%. However, as we will see, for large m K the 
a pole splits in two virtual poles below threshold, one of them moving towards zero 
and approaching the Adler zero region, where the IAM fails. Thus, we will use for our 
calculations the mlAM, although it is only relevant for the mentioned second o pole, 
and only when it is very close to the Adler zero. 



RESULTS 

By changing m n in the amplitudes we see how the poles generated with the IAM 
evolve. We will use the LECs values 10 3 /^ = 0.8 ± 3.8 and 10 3 /£ = 6.2 ± 5.7 from [H 
and fit the mlAM to data up to the resonance region to find 10 3 /[ = —3.7 ±0.2 and 
10 3 /£ = 5.0 ±0.4. These LECs are evaluated at \x = 770 MeV. 

The values of m n considered should fall within the ChPT range of applicability and 
allow for some elastic %% regime below KK threshold. Both criteria are satisfied if 
m% < 500 MeV, since 51/(3) ChPT still works with such kaon masses, and because 
for m n ~ 500 MeV, the kaon mass becomes ~ 600, leaving 200 MeV of elastic region. 

Fig. Q] (left) shows the evolution of the o and p pole positions as m n is increased. 
In order to see the pole movements relative to the two pion threshold, which is also 
increasing, all quantities are given in units of m n , so the threshold is fixed at y/s = 2. 



Both poles moves closer to threshold and they approach the real axis. The p poles reach 
the real axis as the same time that they cross threshold. One of them jumps into the first 
sheet and stays below threshold in the real axis as a bound state, while its conjugate 
partner remains on the second sheet practically at the very same position as the one in 
the first. In contrast, the o poles go below threshold with a finite imaginary part before 
they meet in the real axis, still on the second sheet, becoming virtual states. As m K is 
increased further, one of the poles moves toward threshold and jumps through the branch 
point to the first sheet and stays in the real axis below threshold, very close to it as m K 
keeps growing. The other o pole moves down in energies further from threshold and 
remains on the second sheet. This analytic structure, with two very asymmetric poles in 
different sheets for a scalar wave, could be a signal of a prominent molecular component 
(Kid]]]. Similar pole movements have been also found within quark models lfl2ll . 

Note that ChPT has ben used to evaluate the subtraction constants at the Adler zero 
and the low energy part of LC, always far from the resonance poles, even when they 
move below threshold, as shwon in Fig. Q] (left), being the use of ChPT perfectly justified. 

In Fig. [2] (left) we show the m K dependence of M p and M a (defined from the pole 
position y/Sp i e = M — zT/2), normalized to their physical values. The bands cover the 
LECs uncertainties. We see that both masses grow with m n , but M a grows faster than 
M p . Below m n ~ 330 MeV we only show one line because the two conjugate o poles 
have the same mass. Above 330 MeV, these two poles lie on the real axis with two 
different masses. The heavier pole goes towards threshold and around m n ~ 465 moves 
into the first sheet. Note also that the m K dependence of M a is much softer than suggested 
in fl 1 3M ~ shown as the dotted line, which in addition does not show the two virtual poles. 

In the right panel of Fig. [2] we show the m K dependence of Y p and Y a normalized to 
their physical values, where we see that both widths become smaller. We compare this 
decrease with the expected reduction from phase space as the resonances approach the 
%% threshold. We find that T p follows very well this expected behavior, which implies 
that the pKK coupling is almost m n independent. In contrast Y a shows a different behav- 
ior from the phase space reduction expectation. This suggest a strong m n dependence of 
the o coupling to two pions, necessarily present for molecular states |llll. ll4|l . 

Fig. ED (right) is a comparison of our results for M p with some lattice results |Q]], 
which deserves several words of caution. Our M p is the "pole mass", which is deep 
in the complex plane, and, due to the momentum discretization induced by the finite 
lattice volume, the minimum energy with which pions are produced is larger than the 
resulting M p , so the lattice rho has no width. We can mimic in our formalism a narrower 
p by increasing the number of colors, N c [15]. We also show the result for the rho mass 
for N c = 10. We see that making the p artificially narrower yields a better agreement 
with lattice data. With these caveats in mind our results are in qualitative agreement 
with those of the lattice. Following (|l6l] one may write M p = M p + c\ + 0(m^), where 
the ci parameters are expected to be of order one and M p ~ 0.65 — 0.80 GeV, which 
is confirmed by a fit to lattice data |l6|] . We can fit our results and predict M p = 

0.735 ±0.0017 MeV and c\ = 0.90 ±0. 17 GeV" 1 . Although the pion mass dependence 
of our calculation is steeper than that of the lattice, the values obtained are still consistent 
with the expectations mentioned above. Let us remind that the m K dependence in our 
approach is correct only up to NLO in ChPT. 




FIGURE 2. m K dependence of resonance masses (left) and widths (right) in units of the physical values. 
In both panels the dark (light) band shows the results for the a (p). The width of the bands reflects the 
uncertainties induced from the uncertainties in the LECs. The dotted line shows the a mass dependence 
estimated in Ref. flUl . The dashed (continuous) line shows the m % dependence of the a (p) width from 
the change of phase space only, assuming a constant coupling of the resonance to nn. 
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